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Abstract 

A method is presented to obtain local unitary invariants for multipartite quantum systems 
^^ ' consisting of fermions or distinguishable particles. The invariants are organized into infinite 

Cn I families, in particular, the generalization to higher dimensional single particle Hilbert spaces is 

straightforward. Many well-known invariants and their generalizations are also included. 

1 Introduction 

(— I , The possibility of entanglement between subsystems is a purely quantum mechanical phenomenon, 

^ i' related to the nonlocal nature of the fundamental laws governing our world. In order to understand 

quantum entanglement, one constructs functions on the space of quantum states invariant under 
C ' groups modelling local operations. The two main approaches consider local unitary (LU) operations 

2 . which describe local transformations that can be applied with probability one, or the SLOCC group 

^yi' (stochastic local operations and classical communication) corresponding to transformations that can 

be done with nonzero probability ^4. 

Various results exist for both groups acting on composite quantum systems with distinguishable 
constituents [2", "3", "T, but much less is known about the entanglement of indistinguishable particles. 
»^ ' Recent work revealed that an understanding of fermionic entanglement can also provide us with 

\^ , information about the entanglement of distinguishable subsystems |51 [3]. Motivated by this, in this 

paper we would like to present a way to construct fermionic entanglement measures for pure states. 
Generally, one looks for invariant functions that arc polynomial in the coefficients and their 
C^^ ' conjugates of the pure state with respect to a fixed orthonormal basis. More abstractly (and 

^^ , independently of basis choices), one would like to find subrcpresentations of S{V (BV*) isomorphic 

to the trivial one or one dimensional subrcpresentations of S(V), where V is a representation of 
some group G (either LU or SLOCC) and S{-) denotes the symmetric algebra on a vector space. 

Here, we take a slightly different approach, and use the projections to every subrepresentation 
in S{V) and associate invariants to them. Here V is the state space of a fc-fermion system, and the 
'^ ■ group considered is the LU group, i.e. the unitary group acting on the one-particle Hilbert space. 

5^ , This approach has the advantage that the resulting formulae are independent of the dimension of 

the one-particle state space. 

The outline of the paper is as follows. In section [2] we recall how each graded part of the 
symmetric algebra of a Hilbert space comes equipped with an invariant inner product. In section [3] 
this inner product is utilized in order to associate LU-invariants to every isotypic subspace of the 
space of degree m polynomials in the coefficients of a multi-fermion state. In section U a special 
case is considered, namely, the invariant associated to the subrepresentation containing the weight 
space with highest weight. The method to obtain explicit formulae is also presented here. In 
section [5] some examples are worked out illustrating various features of our approach. In section |6] 
the relationship between LU-invariants generated this way and SLOCC-invariants is highlighted. 
In section [3 it is briefly mentioned, how the fermionic invariants obtained can be used to construct 
local unitary invariants for quantum systems with distinguishable constituents. For the readers' 
convenience, a summary of some concepts from the representation theory of the unitary groups can 
be found in the Appendix. 
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2 The symmetric algebra of a Hilbert space 

Throughout this section % denotes a finite dimensional complex Hilbert space with inner product 
(•, •). We regard "H as a representation of U{J-L) ~ {f -.V. ^ "HlVw G H : \\fv\\ — \\v\\}. 

Let S{H) denote the symmetric algebra of T-L that is, the algebra of polynomials in vectors of 
v.. S{'H) has the structure of a graded algebra, its degree m homogenous subspace will be denoted 
by S"^{'H). As S^CH) — H, and this subspace generates S'('H) as a unital commutative algebra, we 
have that U{T-L) acts on S{T-l) with algebra automorphisms. 

The inner product on "H induces one on S"^{'H) by the following requirement: for u,v d H 
let {u"\v"^) = {u,v)"\ This turns out to be equivalent to saying that for a unit vector u G H, 
\\u"^\\ = 1. Clearly, this inner product will be preserved by the action of UCH) on S{'H), restricted 
to each homogenous subspace. It is known from the representation theory of the unitary groups 
that in this way each S""('H) becomes an irreducible unitary representation of U{H), and hence the 
induced inner product is essentially the only one invariant under this group action. 

To be more explicit, if we fix an orthonormal basis {ei, . . . , e^} in "H, then S"^{'H) is the space 
of degree m homogenous polynomials in the basis elements, and the degree m monomials with 
coefficient 1 form a basis. These monomials are mutually orthogonal, but they are not unit vectors. 
If w = Y,i=i ^i^i then 
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(where L ^"^ fc ) ^^ ^^^ multinomial coefficient) hence 



fci,...,/Cd>0 
ki-\-...-\-kci—m 
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\\vn\'= E (fc,,fc2,...,fcj l"l'"2^---"^'l'll'^l'^2^---^^'ll' (2) 



Comparing this with 
/ d 
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(ikipr- E^n = E (u , ,M'-'.'---'/\' (3) 



we conclude that 



|et^e^^...e^1| = 



m 
fci,fc2,. . . ,kd 



-1/2 



(4) 



3 Invariants for multi-fermion systems 

In this section, H will be a finite dimensional complex Hilbert space, playing the role of the single- 
particle state space of a fermionic quantum system of k particles. If n = dim'H, then the fc-particle 
Hilbert space is isomorphic to 

k k 

/\n^/\c- (5) 



and hence its dimension is (T) ■ This space also conies equipped with an inner product induced 
from that of "H, and a unitary action of U{'H) which models local unitary transformations of the 
fc-particle states. 

Now let us look at the symmetric algebra of the fc-fcrmion state space. On its homogenous sub- 
spaces S"^{/\ H) we have an action of U{7i) which factors through U{/\ H) and an inner product 
which is invariant under U{/\ %) hence also invariant under U{'H). This time the representation 
of U{T-L) is not irreducible, and S"^{/\ H) can be split into the orthogonal sum of f7(H)-invariant 
subspaces in a non-trivial way: 

k 
5-(/\^) = 0Fa (6) 

A 

where A ranges over the partitions of km, and V\ is the corresponding isotypic component of the 
representation. Interestingly, this decomposition is independent of n (apart from the vanishing of 
the subrepresentations associated to partitions involving more than n parts, but for n > km this 
certainly cannot happen) . This is essentially due to the fact that a degree j symmetric polynomial 
in n variables can be reconstructed even if we only now its restriction to a subspace in which only 
j variables take nonzero values. 

This decomposition allows us to introduce unitary invariants, one for each isotypic subspace. 
Let V' G A "H be a /c-fermion state vector, and V"™ its m-th power which is an element of 5""(/\ H). 
Let P\ : S'^{/\ T-L) -^ Vx denote the orthogonal projection. This commutes with the representation 
of C/(-H), therefore the value of /a(V') := {ip"^ , Pxip"") = |1PaV''"IP is invariant: 

V.9 e U{n) : {{g ■ ^'T , Px{g ■ ^)") - {g ■ (^"),5 • (PaV-™)) = (^"\ ^a^™) (7) 

Note that the number of linearly independent invariants is one less than the number of nonvanishing 
isotypic components, because 



A 



^(v™, pav-") = (^"\ \i^Px\i^n = {r\ r^) = i (s) 



Unfortunately, it is in general not an easy task to calculate the projections for all these invariant 
subspaces for every value of k and m, but some of them are easy enough to be done by hand. 

4 Invariant subspaces with maximal highest weight 

Let us now fix an ordered orthonormal basis (ei, . . . , e„) in "H. This also gives the isomorphisms "H ~ 
C", and U{H) ~ U{n, C). The maximal torus T which acts diagonally in this basis is then identified 
with the subgroup of diagonal unitary matrices. The set of one dimensional representations T ^ C^ 
is a commutative group isomorphic to Z". We will use the following identification: 

(ri,r2,...,r„):r^C^ 

" ^ (9) 

(ri,r2,...,r„)(dm5(Ai,..., A„)) = jj A[' 

i=l 

On the set of n-tuples of integers we have the usual partial ordering: (ri, r2, . . . , r„) is called positive 
iff ri -I- ... -I- rn = and ri, ri + r2, ■ ■ ■ , ri + r2 + ■ ■ ■ + r„_i are nonnegative, and A > /i iff A — /x 
is positive. A finite dimensional representation of U{H), when restricted to T, splits into one 
dimensional subrepresentations. The representations with nonzero multiplicity are called weights, 



and a vector whose orbit under T spans a one dimensional subspace is called a weight vector. The 
isomorphism class of an irreducible representation of U{'H) is determined by its highest weight. 
For / ~ {ii, 12, ■ ■ ■ , ik\ where 1 < ii < «2 < • • • < *fe < 't- let us introduce the following notation: 



1 

T 



e/ = eii A ei2 A . . . A a^ = —= 2J cr(7i")ei^(i) «) e^^ (2) «) . . . «) ei^(fe) (10) 



where Sk is the symmetric group on k elements, and a : s^ ^ {1, —1} denotes the alternating 
representation. The set {e|ij_..._ij.}|l < ii < 12 < ■ ■ ■ < ik "^^ n} forms an orthonormal basis of 
[\ H, and therefore every fc-fermion pure state can be expressed uniquely as a linear combination 
of these vectors: 

V" = ^ V'/e/ where ^ IV^/P = 1 (11) 

/e(H) /e(H) 

(Here we have used the short notation [n] = {1,2, ... ,n} and (y) denotes the set of /c-elcment 
subsets of [n].) For each m E N, the jnth power of i/i is a vector in S"^{/\ H): 

V'" = X! V'/iV'/2 •■•'0/„ 6/16/2 ...e/„ (12) 

/i,...,/™ 

We would like to find a vector in 5™(/\ H) which generates an irreducible f7('H)-representation. 
In general we cannot say much about all the irreducible subrepresentations, but we always have 
one weight vector, e™2 fc: corresponding to the highest weight, which is easily calculated to be 
(m, m, . . . , TO, 0, . . . , 0) with k nonzero entries. We now have that (C^('H)e™2.... k) '■— ^ is irreducible. 
The next step will be to find an orthonormal basis for W. 

Our first goal will be to find a generating set for T4^ as a linear space, then we can orthogonalize 
it to yield an orthonormal basis. To this end, we will use the fact that W is also an irreducible 
representation of GL(n, C) whose action on S™{/\ H) is defined in the same way as that of U{H). 

In order to find a generating set which is easy to handle, we will look for one that is the union of 
orbits under Sn < GL{n, C) (possibly up to a nonzero multiple) which permutes the basis elements 
of 7i. It turns out that we can require also that the generating set consists of weight vectors. We 
will call sets with these properties good: 

Definition. Let S C 5""(/\ H) be a subset, ei, . . . , e„ an orthonormal basis in H and Sn < U{'H) 
the subgroup which permutes these basis elements. The subset S will be called good (with respect 
to this basis) if it has the following two properties: 

1. The subset 

k 

CS:= [j Cw'ZS"\/\n) (13) 

wes 

is fixed under the action of Sn- 

2. If V is an element of S then if we write w as a polynomial in the vectors {e/}^ ('[njN then every 

index i € [n] appears the same number of times in every term. Or equivalently: u is a weight 
vector for the maximal torus fixing the given orthonormal basis. 

We can immediately see that {e™}/p('["l^ is the smallest good subset containing e™2 k- 

To reach every element in W, we will use the fact that GL{n, C) is generated by matrices of the 
form Uij (s) = id + sEij where Eij is a matrix with a 1 at the intersection of the ith row and the jth 



column, and zeros everywhere else. We need to know how these matrices act on the basis elements 
of /\ ("H). One can calculate using equation ([T0| that 

e/ i il 

hM • e/ = <{ ei + (-l)l^n(^'^")lse/u{.}\b} jel.iil (14) 

The first and last cases are not interesting, but the second one allows us to build our generating 
set step by step starting from the above mentioned elements. Keeping track of the appearing sign 
could cause some difficulty, but we can overcome this by letting eabc... — —ebac... etc. and simply 
substituting j with i without reordering the indices. 

Observe, that when Uij{s) acts on a degree m polynomial in the e/-s, then we get a polynomial in 
s with coefficients in S'™(/\ H). Since W contains this polynomial for any s € C, and it is a linear 
subspace, W must also contain the coefficient of s' for each < / < ttt. (because the non- vanishing of 
a Vandermonde determinant). Using this method, one can calculate in a few steps a generating set 
for the isotypic (in fact, irreducible) subspace corresponding to the highest weight. The following 
lemma shows which terms should one concentrate on: 

Lemma. Let W < 5'™(/\ H) be an invariant subspace and S CW a good subset 

Suppose that w G S and i ^ j are indices such that i does not appear in w when written in the 
monomial basis as above. Then 

a) The coefficients of every power of s in Uij (s) ■ w as a polynomial in s are weight vectors. 

b) If the degree of this polynomial is d then the one dimensional subspaces spanned by the coefficients 
of s^' and s'*^'' are in the same Sn-orbit. 

c) The coefficient of the constant and the leading terms is contained in <CS. 

d) IfCw ~ Ctt-w' for some tt € S'„, then the minimal good subsets containing S and each coefficient 
in the polynomial Uij{s) ■ w or w^-i(jw-i(j)(s) • w' generate the same subspace. 

Proof, a) Uij{s) ■ ei-^ej^ . . . ej^ = (e/^ + sei'){ei^ + se/^) . . . (e/^ + se/^) where I[ is obtained from 
/; by replacing j with i if /; contains j and e/' = else. The coefiicient of s' contains exactly 
those terms in the expansion in which the number of replaced j-indices is I. 

b) d is the (common) number of occurrencies of the index j in each term of w. The coefficient of 
s'^^'^ term is therefore proportional to the image of the coefficient of s*" under the transposition 
swapping ei and Cj. 

c) The constant term is w. 

d) Let TT <^ Sn <U{'H)he an element such that Cw = Ctt • w' . Then 

Cuy(s) • w = Uij{s)Cw 

= Uj-i(s)CTT ■ w' 
= 'LUij[S)7T ■ W 

= C7ru^-i(i)jr-i(j)(s) -w' 

D 

Corollary. If 5 C S'™(/\ H) is a good subset and w £ S such that in each term of w the index j 
appears exactly once and w does not contain the index i, then Uij{s) ■ w G (S) for all s G C. 



Proof. In this case, Uij{s) • w is a degree 1 polynomial in s, therefore, by the lemma above, both of 
its terms are in CS*, hence their sum is in {S). D 

To sum up, we begin with the vector e^ j,, then act on it and the distinct types of obtained 
coefficients of s successively with the matrices Uij{s), as long as we get new types of vectors. Finally, 
we take union of the ^n-orbits of the vectors we have met. This will result in a generating set of 
W. 

Once we have a generating set, we orthogonalize it, and for each vector w in the orthogonal 
set we calculate the value of |(u;,V''")P 11^11"^. Finally, the sum of these numbers is the value 
of our invariant evaluated on the state V'- Explicitely, suppose, that ip = '^jipi^ij ^-nd w = 
X^fc k k Pki,k2,---,kd^i^ ■ ■ ■ 6/''j whcrc d = (^) and Ii, . . . ,Id are the possible fc-element subsets of 
[n], and ki, . . . ,kd run over nonnegative integers such that their sum equals m. Then by equation Q 



(16) 



ki,...,kd 



We would like to remark that if we are to use these invariants as measures of entanglement, 
then, taking into account the constraint (|S]) and the fact that the mth power of a decomposable 
state is always in the irreducible subspace generated by e^ j,, we should use 1 — (V''",-PvkV''") = 
{ip"^,P^^±tp), or the invariants associated to the subspaces other than W. 

If we wanted to calculate the projectors of the other isotypic subspaces, then we simply needed 
to take the orthogonal complement of W, and find the weight vectors corresponding to the highest 
weight, and proceed with it the same way as we did with e^^ f.. 

5 Examples 

5.1 k = m = 2 case 

The first nontrivial case is the space of quadratic polynomials in vectors of the space of two fermions. 
As we have seen, a weight vector with maximal weight is 6^2, therefore W := {GL{n, C)ef2) contains 
e? for 1 < i < j < n. In the next step we let Ufcj(s) act on an element: 

Wfcj('S)(ey) = (cjj + se^kf = e|- + 2sejjejfc + s'^ef^ (17) 

This shows that wc must add CijCik for each triple i,j, k, where the appearing indices are distinct. 
Now as 

uu{s){etjeik) = (gy + seij){eik + seik) = e^eifc + s{eijeik + etjCik) + eijCik (18) 

we also have to add CijCik + e-ikeij for each combination of indices. 

By the corollary after the lemma, we are ready, but it is instructive to verify the dimension of the 
generated subspace. Clearly, {e^j}i<i<j<n U {eijeik\i<i<n,i=jLj<k^i consists of pairwise orthogonal 
elements. The third type in the generating set is {cijCik + e-ikCij} which is seen to generate a two 
dimensional space for each set of four indices, and these subspaces are pairwise orthogonal and also 
orthogonal to the other elements. Therefore, 

dimW^^f +n ^ +2( =^- i (19) 



2 \ 2 V4/ 12 



which is exactly the dimension of the irreducible representation of GL(n, C) corresponding to the 
partition (2,2). 

Orthogonalization needs to be performed only within the two dimensional subspaces, and this 
leads to the vectors CijCki + eikCji and 6.^6/^ + 2eiiejk — SikS-ji ior 1 < i < j < k < I < n. The 
expression for the invariant corresponding to W is therefore (using equation (|16p ) 

-f(2,2)(V') = {lp,Pwtp) 

n 

= E i^'.i' + E E 2i^.,^.fe 

j=jti=jtk 



l<i<j<n 1=1 l<J<fe<n (20) 



In this case, we can also show that W-^ is irreducible. To this end, let us recall that for n = 4 
there exists a degree two S'L(4, C)-invariant over /\ C*, namely, the polynomial in the Pliicker 
relation which is known to be a sufficient and necessary condition of separability. The subrepre- 
sentation generated by this polynomial is the representation indexed by the partition (1,1,1,1), 
therefore this one must appear also in the n 7^ 4 case. As the dimension of this is (2) , and 

dimT^+ Q = !!(!i^il(^^Zi±^ .dim^^(AC") (21) 

therefore W'^ is irreducible, and the unitary invariant associated to it gives a generalization of the 
Pliicker relation. The explicit formula turns out to be simpler than the previous one: 

2 

^(1,1,1,1) (V") = ('0,-Pw-lV') = E -^li'ij'ipki + AkiPij + Aii'jkl'^ (22) 

l<i<j<k<l<n 

5.2 k = 2, m = 3 case 

In this case, a weight vector for the highest weight is 6^2- Again, W :— {GL{n, C)ei2). We are 
looking for a generating set of W. We extend 6^2 into a good set {eij}i<i<j<„. Now we need to 
add the coefficient of s in 

W32(s)(ei2) = (ei2 + seia)^ = e^ + Sscijeis + is'^e^el^ + s^e^g (23) 

and one vector from each element of the orbit of the subspace generated by it: {ef,eik}ij^k£[n]- The 
next steps are: 

U43(s)(ei2ei3) = ei2(ei3 + se^)^ = . . . + 23612613614 + s^(. . .) (24) 



Wml(s)(6i26i3) = (612 + S6m2)(ei3 + SCms) 

(25) 



+ s(26i26i3e,„3 + 6,„2e?3) + s^(. . .) + s3(. . .) 



Here m = 2 is special, in this case the second term in the coefficient of s vanishes, hence we have 
to add {^eijeikGmk + ^mj^ik} for any ordered pair of disjoint pairs (i, j, k, m), and also eijCjkekj for 
{i,j,k} e ( 3 )• The remaining steps are 

Mml(s)ei2ei3ei4 = (ei2 + S6,„2)(ei3 + Se„i3)(6i4 + S6m4) , , 

2 3 ' ") 

= . . . + S(6m2ei36l4 + 6126^3614 + 6126136^4) + S (...)+ S (...) 



Mml(s)(e52ei3ei4 + 612653614 + 612613654) = . . . + 5(6526^3614 + 6526136^4 

+ 6^2653614 + 6126536^4 (27) 

+ 6m26l3654 + 6126^3654) + S (...) 

Here m < 5 does not lead to a new subspace. 

It turns out that the vectors obtained so far are enough to generate W. In this case, orthogonal- 
ization turns out to be a bit lengthy, especially in the case of the six-term vectors like in equation 
((27|) . These span a five dimensional subspace for each six-element set of indices ii, . . . ,iQ. For these 
the coefhcients of the monomials are given as a matrix: 
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-1 
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-1 
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-1 


-1 
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-1 


-1 





-1 


-1 





\4 
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-2 
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-1 


-2 


-1 
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-1 


1 


2 


1 


-1 


-2/ 



(28) 



The order of the monomials is (12)(34)(56), (12)(35)(46), (12)(36)(45), (13)(24)(56), (13)(25)(46), 
(13)(26)(45), (14)(23)(56), (14)(25)(36), (14)(26)(35), (15)(23)(46), (15)(24)(36), (15)(26)(34), 
(16)(23)(45), (16)(24)(35), (16)(25)(34), where {ab){cd){ef) is a short notation for 6,„,,,e^^,i^e,^,i^. 
The norms inverse squared of these vectors are 

13 3 1 

^'8' 8' 8' 8 ^'^^' 

respectively. The orthogonal generators coming from the remaining vectors are given in table ([T]) 



form 


indices 


dimension 


II • ir^ 




{*,J} 


(2) 


1 


6i,-6,fc 


«,{j},W 


n{n - l){n - 2) = QQ 


3 


Sijeikeu 


W,{j,fc,^} 


"(V)=4(:) 


6 


^ij^ik^kj 


{h3,k} 


I3J 


6 


ef-jSki + 2eijeiiekj 


{^.J}AKl} 


2(2) ("2^) 


1 


-2ef^eki + GeijCikeij + 2eijeiiekj 


i 

8 


(^ikf^il^jm ~r G-ik^jl^im ~r (ijk^il^im 


{i}:{j>fc:':"l} 


3-rr) = i5(;:) 


2 


eikeuejm + •ieikeijCmi + etmeackj + 

^^im^ij^kl ~r ^G-imf^ik^-jl 


i 
4 


2eijeiiemk + eikeuemj + eikeijCmi + 

^im^ilGjk ~r (^ivn^ij^kl 


3 
4 



Table 1: Orthogonalized generators for W . Indices shown in one set are indistinguishable for 
counting purposes. 

Using these data, the value of the invariant /(3.3) can be calculated in a straightforward way, 
but the full formula is too long to be presented explicitely. 

The orthogonal complement of W clearly has a highest weight of (2, 2, 1, 1), and we could find 
a generator of the unique one dimensional weight space corresponding to it, and calculate the 
projector of its invariant subspace. Instead of this, we follow another approach. According to the 
plethysm 



3(3)[S(1,1)] = S(3^3) -I- S(2,2,l,l) + 5(1,1,1,1,1,1) 



(30) 



for n ~ 6, an SL{6, C)-invariant polynomial appears. It is easy to guess how this should look like: 
for a state ip E /\ C^, we can construct ip f\ip /\iIj which is an element of /\ C^, a one dimensional 
vector space on which GL{6, C) acts by multiplication with the determinant. Therefore this element 
remains unchanged under SL{6, C), and its norm squared is an f/(6, C)-invariant polynomial in the 
coefRcients of "0 and their conjugates. Our invariant corresponding to the subrepresentation indexed 
by the partition (1, 1, 1, 1, 1, 1) must be proportional to it. Explicitely, it equals to 

2 

1 



11520 



TrGSe 



O'(7r)'07r(l),7r(2)V'7r(3),7r(4)'07T(5),7r(6) 



(31) 



Here, the sum is over all the permutations, but actually there are 15 different terms, each counted 
48 = 3! • 2'^ times. Alternatively, we could sum over the partitions of [n] into three two-element sets. 
The n > 6 case can be obtained similarly to the previous section. Taking all the six-element 
subsets of [n] polynomials like this span an (") dimensional subspace which is also the dimension 
of the invariant subspace we are looking for. Therefore in the general case the invariant is 

2 



'(1,1,1,1,1,1 



m = 



11520 '^ 

/e(H) 



E 



'^i''^)4'i„{l),i„i2)^i„{3)A„{4)'4'i^(5),i„{6) 



(32) 



where / = {ii, . . . ,«6}- 

These two invariants are linearly independent, and they sum to 1 with the one associated to the 
third irreducible subspace. 

5.3 k = 3, m = 2 case 

Now we turn to the first case with more than two particles. In S^{/\ H) the vector with highest 
weight is 6^23- We proceed in a similar way as before: 

Mn3(s)(ei23) = 6^23 + 2sei23ei2n + s^e^jn (33) 



W„2(s)(ei23ei24) 



+ s(ei23ei„4 + ei„3ei24) + s^(- ■ •) 



(34) 



Wnl(s)(ei23ei54 + 61536124) = ■ ■ . + s(ei236„54 + e„23ei54 + ei536„24 + e„53ei24) + s^(. . .) (35) 

These vectors already form a generating set, we only need to orthogonalize this set. For a fixed 
subset of six indices, the vectors of the form like in (j35p span a five dimensional subspace. Orthogonal 
generators for this are again given with the coefficients of the monomials as a matrix: 



/ 1 



V 





2 



2 

-1 



1 
1 

-1 
-1 
-1 





2 
2 
-1 





2 


-2 
1 





2 

-2 
1 



1 \ 
-1 
1 
-1 

-1/ 



(36) 



The order of monomials is (123)(456), (124)(356), (125)(346), (126)(345), (134)(256), (135)(246), 
(136)(245), (145)(236), (146)(235), (156)(234), where {abc){def) is a shorthand notation for the 



(37) 



vector ei^^i^,i^ei^^i^^ij. The inverse squared norms of these vectors are 

11111 

2'6'6' 18' 9 
respectively. The orthogonal generators coming from the remaining vectors are given in table ^ 

The value of /(2,2,2) can now be calculated. This time the orthocomplement is also irreducible, 
so we get one independent invariant in this case. 



form 


indices 


dimension 


■ w^"^ 




{»,J,fc} 


(3) 


1 


^ijk^ijl 


{^.J}Ak,l} 


(")(V) 


2 


^ijm^ikl 1 f^ijk^iml 


{i}-,{hk,l,m} 


•^^T) 


1 


^ijm^ikl 1 ^^ijl^imk ^ijk^iral 


i 
3 



Table 2: Orthogonalized generators for the subspace generated by the highest weight vector. Indices 
shown in one set are indistinguishable for counting purposes. 

6 SLOCC-invariants and local unitary invariants 

In the examples we have seen local unitary invariants with a special property: for a particular 
value of n, the corresponding irreducible subspace becomes one dimensional, and the subspace is 
pointwise fixed under the action of SL{n, C), that is, the local unitary invariant turns out to be a 
SLOCC-invariant. Let us examine this case in more detail. 

The irreducible polynomial representation of SL{n, C) indexed by the partition A is one dimen- 
sional precisely when A consists of equal parts. In this case, A = (r, r, . . . , r) is a partition of nr, 
hence a neccesary condition for it to occur as a subrepresentation of S'™(/\ C") is that mk = nr, 
and in this case GL{n, C) acts on it by multiplication with the rth power of the determinant. The 
norm squared is therefore invariant under U{n,C). 

In our notations, this subspace is spanned by a polynomial w in the basis vectors ei, . . . ,e„. 
w is a weight vector with weight {r,r, . . . ,r), and it generates a one dimensional U{n, C)-invariant 
subspace. The crucial thing is that when we increase the dimension n of the single particle state 
space to n' , w remains a weight vector that generates an irreducible C/(n', C)-invariant subspace, 
but it is no longer one dimensional. Therefore, the invariant corresponding to this subspace will be 
a generalization of the SLOCC-invarant we have begun with, but is now only a unitary invariant. 

The explicit form of the resulting invariant can be obtained in general using the method outlined 
above: we must act on it with Uij(s)-s and elements of Sn'- A particularly simple special case is 
when r — 1. In this case the dimension of the representation corresponding to A is ("j), and an 
orthonormal basis can be obtained by acting on w by elements of Sn' ■ Therefore the invariant can 
be obtained by calculating the value of the SLOCC-invariant with the initial index set [n] replaced 
by every element of ( " ), and summing their absolute values squared. 

7 Distinguishable particles 

We can also obtain many (but not all) local unitary invariants for distinguishable subsystems from 
our fermionic ones [71 [3 [6]. Suppose that we would like to find local unitary invariants for a quantum 
system containing k subsystems with Hilbert space dimensions ni, . . . , n^, and let n = ni + . . . + nk. 
Then using the branching rule 



U{n,C) > U{ni,C) x ... x U{nk,C) 

k k 

/\(n)^/\((ni)©...©(nk))^.. 



(ni,n2,...,nk) 



(38) 



we can identify the full state space with a subspace of a fermionic Hilbert space with k particles 
and n dimensional single particle state space. Then we can pull back the fermionic invariants 
constructed by the method above. For example, for fc = 2, the restriction of /(i, 1,1,1) is proportional 
to the square of the norm of the concurrence vector [3]. 

Observe, that these special LU-invariants have a larger symmetry group, in particular, if ni = 
. . . = rife, then they are also permutation invariant. 
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As an other example, let us consider the quantum system of three qubits. Using the scheme 
outlined above, its Hilbert space can be thought of as a subspace of the Hilbert space of a three- 
fermion system with six single particle states. In the previous section (see equation ([36]) and table[2|) 
we have derived a formula for a local unitary invariant which can be restricted to this subspace. 
The restriction turns out to be the following: 

1 1 

^(2,2.2) W = Yl I'^y-fcl^ + 2 X! (l^uoV'jilP + li'iOjAljl'^ + IV'OyV'Hjn 
t j',/c— * J— 

1 / I 

+ y^ IV'iOlV'ilOp + ol^«l^jlO + 2V'jOoV'illP + IV'OilV'HOp 

+ IV'OilV'HO + 2V'OioV'lilP + IV'OliV'lOiP + o IV'OljV^lOi + 2tp00i1plli\'^ ] (39) 

+ ol^oooV'iiiP + c IV'oooV'iii + 2i/'ooiV'iioP + ^|'0oooV'iii + 2i/'oiiV'iooP 

2 D D 

+ T^IV'oooV'iii - 2V'ooiV'iio - 2^oiiV'iooP 

+ q IV'oooV'iii + V'ooiV'iio + Si/'oioV'ioi + V'oiiV'iool^ 

where ip = ^^ ■ ^^q ipijkSi ® Cj ® e^ is a three-qubit state. This quantity is invariant under local 
unitary transformations and permutations of the three subsystems. Plugging in the coefficients of 
certain states into this formula reveals that for a separable state the value is 1 as expected, for 
the GHZ-state it is |, for the W-state |, and for a biseparable state with maximal entanglement 
shared by two qubits, the value is |. This behaviour is in contrast with the well-known three-qubit 
invariant, Cayley's hyperdeterminant, which is only sensitive to GHZ-type entanglement. 

Note that these values are one minus the quarter of the values of the entanglement measure 
proposed in [^ for three qubits, so one is tempted to conjecture that the restriction of 4 — 4/(2,2.2) 
equals to their measure. In their paper this measure appears as a member of an infinite family of 
multiqubit entanglement measures. It would be interesting to see whether every member can be 
found using our method. If this is the case, then we could generalize them to arbitrary dimensional 
single particle states and also to fermionic systems. 

A more detailed treatment of this method of constructing entanglement measures for distin- 
guishable subsystems from fermionic measures along with more examples can be found in [B]. 

8 Conclusion 

In this paper a certain class of entanglement measures for fermionic quantum systems has been 
introduced and studied. A way to obtain their explicit form is presented, and it was pointed out 
that this form is independent of the dimension of the single particle state space. Some examples are 
discussed in detail. The connection to SLOCC-invariants and the case of distinguishable subsystems 
is also mentioned. 

At this point a number of questions arise. Further study is needed to explore the behaviour of 
these local unitary invariants. For instance, are there any entanglement monotones among them, 
and is there a way to characterize these? Could we find the convex roof extension of any of them? 
Is there a physical meaning for these quantities, in what way do they measure entanglement? 

Also, similar method could be applied directly to quantum systems with distinguishable parti- 
cles, or to mixed states. It would be interesting to see if one could obtain entanglement measures 
this way which can be useful in practice. 
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9 Appendix: Representation theory of the unitary group 

In this appendix, the basic aspects of the representation theory of the unitary groups is summarized. 
These and many more can be found in many textbooks, see e.g. [9]. 

The unitary group C/(n, C) consists of n x n complex matrices satisfying A* — A^^. Let 
T < U{n,C) be the subgroup of diagonal matrices, T is a maximal torus in C/(n, C). We are 
interested solely in finite dimensional polynomial representations of U{n, C) which are determined 
by their characters which are in turn uniquely encoded in the restriction of the characters to T. 
The characters themselves are symmetric polynomials in the eigenvalues (the diagonal elements, in 
the case of T) . 

For example, the action of U{n, C) on the n dimensional vector space C" of column vectors by 
left multiplication is called the standard representation. Its character is represented by symmetric 
polynomial xi + X2 + . . . + x„. 

A representation V of U{n, C) is also a representation of T whose irreducible representations 
are one dimensional, and hence can be considered to be homomorphisms T — >■ C^. An irreducible 
representation of T with nonzero multiplicity in V is called a weight. The weight space corresponding 
to weight is the union of the subrepresentations in V isomorphic to a given weight. 

Continuing the previous example, the weights of the standard representation are the represen- 
tations pk ■ diag{Xi, . . . , A„) i— ?> A^, and the weight space of pk is spanned by the fcth standard basis 
element. From now on, the product p^^ (g) . . . (g) pl^" will be simply denoted by (ri, . . . , r„). 

There is the usual partial ordering on the set of weights which we identify with Z": (ri, r2, . . . , r„) 
is called positive iff ri + . . . + r„ = and ri, ri + r2, ■ ■ ■ , ri + 7-2 + . . . + r„-i are nonnegative, 
and A > /x iff A — /i is positive. The isomorphism class of an irreducible representation of U{n, C) 
is determined by its highest weight and the weight space for the highest weight is one dimensional. 
Using this fact, we can decompose any finite dimensional representation of U{n,C) into the direct 
sum of irreducible representations. 

The dimension of an irreducible representation of U{n, C) with highest weight A = (Ai, . . . , A„) 
is 

dim^(,,,...,,„)F= n "''-^''^ (40) 

l<z<j<n 

Note that if we start to increase n and pad A with zeros on the right, then the dimension in 
the function of n turns out to be a polynomial of degree Ai + . . . + A„ . In fact the entries of the 
representing matrices are polynomials in the entries of the represented matrix with this same degree. 
The symmetric polynomial giving the character of this representation is the Schur-polynomial. 

Calculating symmetric and exterior powers, tensor products, decomposition to irreducible sub- 
representations and many more can be done working only with symmetric polynomials. This is a 
crucial fact for our purposes, as a symmetric polynomial of degree d is determined by its terms 
containing unknowns only from a fixed set of d variables. Indeed, no term can contain more than d 
variables, so the missing ones are obtained by permutations of the variables. This fact implies that 
the decomposition of S"''{/\ C") to irreducible C/(n, C)-representations is independent of n when 
n < km. 
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